This study develops a numerical steady-state solutions to both buoyancy driven and surface tension driven flow problems in two spatial dimensions. A high-order spectral method is ap- plied for the spatial discretization, ...